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Abstract 


The high-charge high-energy (HZE) ion components of the galactic cosmic rays present unique challenges 
to biological systems in comparison to terrestrial forms of radiations. We present a deoxyribonucleic acid (DNA) 
breakage model to visualize and analyze the impact of chromatin domains and DNA loops on clustering of DNA 
damage from X-rays, protons, and HZE ions. Our model of DNA breakage is based on a stochastic process of DNA 
double-strand break (DSB) formation that includes the amorphous model of the radiation track and a polymer model 
of DNA packed in the cell nucleus. Our model is a Monte-Carlo simulation based on a randomly located DSB cluster 
formalism that accommodates high- and low-linear energy transfer radiations. We demonstrate that HZE ions have 
a strong impact on DSB clustering, both along the chromosome length and in the nucleus volume. The effects of 
chromosomal domains and DNA loops on the DSB fragment-size distribution and the spatial distribution of DSB 
in the nucleus were studied. We compare our model predictions with the spatial distribution of DSB obtained 
from experiments. The implications of our model predictions for radiation protection will be discussed. 


1. Introduction 


O NE of the main concerns for long-term deep space missions is health risks associated with prolonged exposure 
to proton and high-charge high-energy (HZE) components of the galactic cosmic rays (GCR) and solar 
radiation in space. It is known that an astronaut on a six-month journey to Mars, which is the time required with 
conventional propulsion, would be exposed to about 0.3 Sieverts (1 Sievert = 1 Gray = 1 Joule/kg = 100 rad = 100 
rems for X rays, but is equal to Q X 1 Gray =100 rems = QxlOO rads for high-linear energy transfer (LET) radiation, 
where Q is the biological quality of radiation; Q > 1 for high-LET), or even to 0.6 Sieverts on a round-trip 
(Cucinotta et al., 2005). Eighteen months on the surface (if it takes that long to get there, you might as well stay for 
a while) would bring another 0.4 Sieverts, for a total exposure of 1 Sievert. Limits set by NASA vary with age and 
gender ranging from 0.5 to 2 Sieverts.. During a 3-year mission to Mars at solar minimum (worst case for the GCR 
exposure), 46% of cells might be hit by an HZE particle (with the electric charge Z > +15), with 13% hit by an iron 
particle (Z = +26). Hit frequencies are 4-8 times lower for nuclei of cells. The cell nucleus would also be traversed 
by a proton twice a week, and an a-particle once a month (Badhwar et al., 2002; NCRP report 98; Townsend et al., 
1992; Cucinotta et al., 1998b). Particle fluences from GCR may be more relevant than absorbed radiation doses, 
which will be a few tens of cGy, to the human body. This is because the radiation effects measured at short times 
after fractionated doses of 100 cGy or higher, 
which correspond to several hits per cell, may 
not accurately predict astronaut risk because of 
deoxyribonucleic acid (DNA) repair. Repair 
and adaptation mechanisms (Rothkamm & 

Lobrich, 2003), which can counteract effects of 
particles delivered chronically, may be over- 
whelmed by the delivery of radiation as an 
acute dose from one fraction. It is also impor- 
tant to remember that other conditions during 
spaceflight may modify responses to HZE 
radiation. 

GCRs are the major radiation source in 
outer space, and these are composed of H (85%), 

He (13%), and HZE particles (1%). The relative 
abundance and ionizing power of the even- 
numbered HZE particles in GCR is well studied 
(Report LBNL-40278). Although the fluence of 
HZE particles in GCR is low, the high charge of 

these particles makes them particularly damaging (Fig. 1). Because of its relative abundance, and charge of Z = +26, 

Fe is a major concern among the various HZE particles encountered. Although C, O, B, Li, N, and Be particles are less 
highly charged, their abundance indicates they must also be considered. 

Some knowledge about exposure to HZE particles in outer space has already been gained, as the HZE dose 
rate during the Apollo 17 flight was calculated to be 0.12 pGy per hour, or about 1.05 mGy per year for HZE particle 
irradiation. Ignoring solar particle events (SPEs), the total dose of radiation during a 3-year mission is about 0.3 Gy, 
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Charge Number 


Figure 1 . GCR charge contribution. 
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mostly due to protons. It has been estimated that during a Mars mission, annual exposure from GCR would be between 
0.15 Sv/year at solar maximum and 0.58 Sv/year at solar minimum. It has been estimated that SPE exposure (from 
electrons and protons as well as other secondary radiation induced by them) may contribute about another 0.3 Sv/year, 
or almost 1 Sv on a 3-year Mars mission, although the exposure rate due to an SPE could increase 100- to 1000-fold 
during brief intervals (up to several days). However, in contrast to GCR, shielding can effectively mitigate SPE doses 
with their energies and ranges in materials being much lower compared to GCR. 

Energy deposition from GCR is largely confined to a thin cylinder of tissue that receives a high local dose, 
especially at the end of the particle range, which is the depth of tissue penetration (Turner et al., 2002). One can 
calculate that for a 1-GeV/u 56 Fe particle, the average dose to a nucleus that would be traversed near the nucleus center, 
as recent data by our model show, would be about 12 cGy. Furthermore, the effect of ionizing radiation on cells and, 
ultimately, tissues originates primarily from damage to DNA; and the complex ways in which particular types of ra- 
diation interact with matter and break and/or otherwise alter DNA structures control the potential consequences of 
radiation, modulated by the ability of cells to repair damage. The passage of an HZE particle through the cell nucleus 
should cause multiple, intense, and essentially instantaneous ionization events that induce complex patterns of DNA 
damage, sometimes irreparable. Little is known about the effects of charged particles at the cellular and molecular level 
in mitotic cells. Mitochondrial as well as nuclear DNA can be a radiobiological target. Since the mitochondrial electron 
transport chain is the main endogenous source of reactive oxygen species that causes oxidative stress, damage to mito- 
chondrial DNA can be particularly relevant during spaceflight. Mitochondria have multiple copies of their genome and 
even possess DNA repair enzymes. However, the studies of the effects of HZE radiation in this area are insufficient. 

The radiation damage of DNA is crucial to understanding the effect of radiation on an astronaut’s body 
(Goodhead, 1995). DNA is a set of long molecules comprising 46 human chromosomes in a normal diploid genome. 
The female genome is slightly longer because of the double X chromosome, which is much larger than the Y chromo- 
some. We have developed computer models, and related codes, as well as visualization software of damage to DNA 
from charged particle tracks. For the current analysis we simulated the male genome with XY chromosomes. Our code 
can easily allow for the double X chromosome. The genomic size is typically measured in the number of base pairs 
(bps) of the well-known double-helix molecule, or its derivatives, kilo base pairs (kbps) and mega base pairs (Mbps). 
The genomic size is proportional to the molecular weight of the chromosomes, even though it far exceeds the molecular 
weight of the bases, containing the genomic information, because chromatin, or the fiber that the chromosomes are 
made of, is a complex collection of structural proteins, mostly (Rogakou et al., 1998) histones (Fig. 2). The total length, 
or molecular weight, of the human genome, is about 3.2 Gbps (giga base pairs). 

Radiation causes ionizations in cell water that leads to aggressive water radicals attacking DNA. One consequence 
of radiation is the loss of bases that have the genomic information (Taucher-Scholz, 2000). The loss of one base 
is called a single strand break. If two single-strand breaks occur within about 10-50 bp, a DSB is produced (Wu et al., 
1996; Brenner & Ward, 1992). Although 
methods for measuring DNA damage have been 
developed over the years, it is still unclear how 
well damage can be measured at low doses 
because of the small regions or molecular sizes 
of the damaged regions of interests. The most 
common measurement is the pulsed-filed gel 
electrophoresis experiment (PFGE), where after 
special preparation of DNA and the use of an 
applied electric field, a collection of fragments 
produced by DSB is sorted by the molecular size 
and the DNA fragment distribution is measured. 

The number of fragments is almost precisely 
equal to the number of DSBs. The limitation of 
this experiment is that small fragments, as well 
as very large ones, cannot be detected. There are 
very few large fragments, so this number is not 
significant; but the number of small fragments is 
large and dependent on the radiation quality, 
and, thus, is an important consideration in 

determining the effects of space radiation. Here we report our novel method of accounting for missing fragment and 
leading to new estimates of DSB counts induced by radiations relevant to space environment. We also compare our 
approach to other approaches (Results section). 
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Figure 2. The high-order structure of DNA is shown 
on the level of chromatin fiber. The source for this is 
http://sqi.bls.umkc.edu/waterborq/chromat/chroma09.html . 

A piece of DNA depicted is about the size of one monomer that 
will be a building block for the whole genome in this work. 
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Another important consideration is the geometric statistics of DSB (Prise, 1997; Pinto et al., 2004). The same 
total number of fragments may have a different relative number of smaller fragments. PFGE experiments report a num- 
ber of smaller fragments in excess than one would expect if DSB were located randomly (Campa et al, 2004) along 
chromosomes. This is believed to be the signature of the high-LET radiation, as well as the result of the particular 
properties of DNA organization in the cell nucleus. This phenomenon, called DSB clustering, is herein addressed by 
our model. The importance of these data is that DSB clustering indicates a higher severity of damage. A collection of 
closely located small fragments will more likely result in a permanent genomic damage, such as gross violation of 
DNA topology called chromosomal aberrations, and, as related to them, the loss of genomic information (gene 
deletions). 

DNA damage can lead to cell death, which is not problematic for small doses during space missions. (In the major- 
ity of the brain tissues, however, cells do not replenish, so cell death is potentially dangerous.) More importantly, the 
accumulation of genomic damage can lead to late effects, such as cancer, through genomic instability and mutations. 

At the NASA Johnson Space Center, we have developed a suite of numerical and mathematical tools to address 
the biophysical modeling of radiation, cells, and tissue. We will present here a brief description of methods used in 
the present work. First and foremost, we have developed precise physical models of radiation transfer through matter 
(Cucinotta et al. 1998a; Cucinotta et al. 1998b; Cucinotta et al., 1995; Cucinotta et al., 1999; Cucinotta et al., 2000; 
Nikjoo et al., 2002a; Nikjoo et al, 2002b; Wilson et al, 1987). These models address several types of radiations: low- 
LET, such as X rays, and gamma-rays; and high-LET radiation, which are primarily a-particles and heavier in mass, A, 
higher in charge, Z, and higher in energy E, FIZE ions. LET is a significant parameter indicating the quality of radia- 
tion. Its meaning is simple; it is the energy in keV imparted to the matter (such as cell cytoplasm) per pm of the particle 
path. The more energy imparted per pm, the more biologically destructive the radiation is (for very high LET, there is a 
drop-off of the radiation efficiency due to overkill effect (NCRP report 98)). In the present analysis, we treat low-LET 
radiation as a homogenous field of radiation that does not have sharp peaks of deposited energy in small volumes. By 
high-LET radiation, we mean ions such as 1 GeV/a 56 Fe (an HZE ion), or a-particles, that have a sharp peak of deposit- 
ed energy close to the track center. Some researchers argue that HZE ions are not strictly high-LET radiation due to 
their extensive penumbra (energy field at distances up to 50 pm from the track core), where energy is imparted by 
low-LET electrons (5-rays). However, the track core (a narrow region around the particle trajectory) is high-LET, 
and it is more significant for the biological effects as our studies have shown. 

We should distinguish between LET, which is an average measure of the energy loss along a track in z 
direction and integrated over the span of the track in x and y directions, and the track’s energy profile in the x-, y- 
plane perpendicular to the particle direction. Often higher LET particles such as 56 Fe have a wider penumbra, or the 
extent of the energy field, than the particles with low LET. For instance, a 151-keV/pm 56 Fe ion deposits energy in 
the cylinder of up to 2*30 pm in diameter and more, while 40-keV/pm a-particles have a track confined to less than 
0.1 pm of the high-energy deposition area, almost like a needle piercing cells. This field of energy surrounding im- 
pinging particles is simulated to a high precision by the program that we use in this work (Cucinotta et al, 1999; 
Cucinotta et al., 2000). 

We have addressed the effect of radiation on cellular level, and designed a sophisticated model of DNA 
breakage in a cell (Ponomarev et al., 2000). This model includes the track structure described above and a polymer 
model of chromosomes based on a random walk (RW) model for chromosomes (Fig. 3). This was a significant attempt 
to address the biological effects of radiation insult on cells based on our state-of-the-art knowledge of the underlying 
physical processes. 

One of the main motivations for radiobiological research is to guide extrapolations of risk estimates from 
high to low radiation doses and associated dose-rates. By “low” we mean the dose-rate encountered in space, while 
the high dose-rate is of typical radiobiological experiments that have practical limitations on what dose and dose- 
rates conditions can be used. This is of importance for individuals environmentally exposed (e.g., from radon) or 
occupationally exposed (e.g., air-flight personnel and astronauts) to low- or high-LET radiations. For example, 
recent PFGE experiments (Lobrich et al., 1996; Hoglund et al., 2000; Stenerlow et al., 2003; Belli et al., 2000) provide 
DNA breakage yields at doses that are too high to be directly applicable to radiation risk predictions in humans. 
However, if a reliable model allowing extrapolation to low doses is available, high-dose data can provide the basis for 
useful information about DNA breakage yields at low doses. For low-LET radiation, the random-breakage model 
(Blocher, 1990) is suitable for such a task. 

The random-breakage model (Cook & Mortimer, 1991) that was used ubiquitously in the past is not adequate for 
high-LET radiation because it lacks mechanistic detail to allow for the pronounced clustering of DNA damage caused 
by the highly inhomogeneous energy profile of high-LET ion tracks. Specifically, the DNA is multiply folded during 
interphase; but, even on large scales, DNA loci, which are close to each other along the DNA contour, tend also, on av- 
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erage and neglecting fine- structure, to be spatially close to one another. With this geometry, a bias for ionizations to be 
close to each other in space at high LET produces a bias for DSBs to be close to each other along the DNA. 

A new model, called NASARadiationTracklmage (Ponomarev et al., 2000; Ponomarev et al., 1999; Ponomarev el 
ah, 2001a), based on radiation track structure and large-scale chromatin geometry has been shown to predict the non- 
randomness of high-LET-induced DSB distributions within the genome, in the NASARadiationTracklmage model, ion 
track structure is simulated (Cucinotta et al, 2000) and is superimposed on random-walk (Sachs et al., 1995) chromatin 
geometry (Fig. 3a). To make the picture more realistic, we introduce chromatin loops (Fig. 3b), which are found in ex- 
periments revealing high-order DNA organization (Marco & Siggia, 1997). Ionizations due to a high-LET radiation 
track are spatially correlated, being predominantly near the line representing the center of the track (Fig. 3c) rather than 
spread randomly over a whole cell nucleus. This clustering of ionizations, which is determined by the type of incident 
particles, leads to correlations among DSBs along chromosomes as opposed to the mostly uncorrelated DSB produced 
by low-LET radiations. The advantage of this approach is that, since it is a generalization of the random-breakage 
model (suitable for X -rays, for example), it is easy to reduce the NASARadiationTracklmage model to the case 
with low LET (Ponomarev et al., 200 lb). 



Figure 3a. The NASARadiationTracklmage graphics user 
interface (GUI) displaying one Monte-Carlo realization of RW 
(each chromosome given in a separate color), a track center 
(red line), and a round nucleus (the outline of the membrane is 
intuitively visible). Overlap between RWs may or may not be 
allowed. The software is able to simulate many tracks, an 
ellipsoidal nucleus, RW packed into defined chromosomal 
domains, certain impact parameter for the track, and a fixed 
fluence of particles. Each chromosome depicted here is a vast 
collection of dots (about 100,000) corresponding to monomers. 
Here, the scale of the whole nucleus individual dots is not 
discernable. The size of a chromosome corresponding to 
one dot is, however, about the size of a piece of DNA and 
the other accompanying molecules shown in Fig. 2. 
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Figure 3b. A group of 10 loops is shown on a much smaller 
scale than in Panel A. One chromosome can be divided into 
10-100 such structures, each structure containing 1320 kbps 
of genomic material. Loops in the model are produced by 
random walks that return to the origin. 
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Figure 3c. A group of DSBs is shown, as well as an Fe track, 
and a chromosome #6 (other chromosomes are present but 
are not shown). DSBs are located predominantly along the 
track center. Due to an extended penumbra of the Fe ion 
track, some DSBs can form at a distance from the track center. 
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The proposed mechanistic approach predicts DSB fragment-size patterns in the genome over a very large range 
of fragment sizes (20 kbp-10 Mbp); this range has been successfully compared to the experimental data for a variety of 
situations (Ponomarev et al., 2000; Ponomarev & Sachs, 1999; Sachs et al, 1995; Sachs el al., 1999) where distribu- 
tions of DNA fragment sizes were measured after irradiation at LET=80 keV/pm or more (Lobrich et al, 1996; 
Hoglund, et al., 2000; Newman et al., 1997). 

The NASARadiationTracklmage model was previously formulated for a simple polymer model of DNA that 
had only RW for the chromosomal structure. In this work, we introduced the chromosomal domains and DNA loops 
as a more precise model of the high-order DNA organization in the cell nucleus. Our analysis of the impact of these 
structures on the fragment sizes is presented herein. 


2. Methods 

2.1. The application of randomly located clusters formalism 

Our polymer model of DNA breakage is based on a more general ( i.e ., than the random breakage model) 
randomly located clusters (RLCs) formalism (i.e., the clusters ofDSBs). The generality of the formalism is based on 
the fact that we do not need to know a priori the reason for DSB clustering, or the precise form of it, to infer breakage 
from many particles (high dose, or high fluence). The determination of the unknown function of clustered DSB can be 
delegated to a Monte -Carlo algorithm that we will describe below. Detailed descriptions of the RLC formalism and the 
Monte-Carlo algorithm employed in NASARadiationTracklmage are given in Ponomarev et al, 2000, and in Sachs et 
al, 1999, along with the validation and the benchmark of the algorithm. The fundamental ideas underlying the RLC 
formalism are that the one-track action can make a stochastic cluster of DSB along a chromosome; that different clus- 
ters on a chromosome, due to different tracks, are statistically independent; and that the location of DSB clusters in the 
genome, unlike the non-random location of DSB within one cluster, is random, or Poisson (hence the term RLCs). 

The RLC formalism relates low-dose, one-track effects to high-dose, multi-track effects. Thus, only the one- 
track action needs to be simulated. This holds true for various radiations and various chromatin geometries, as long as 
the spatial distribution of incident ions traversing a certain area is random. It is assumed that there is no change of LET 
along the ion trajectoiy on the scale of a cell nucleus. The main result of the formalism is the inference of the multi- 
track action through the one-track action, as given by the following equation, stated for brevity for a genome with 
chromosomes of equal size S (the generalization to a more realistic genome is straightforward and is done 
computationally) : 


1 ~ R(s) - [l - s/5]/ (s) exp 


-A 


s + 


(Af-l)j(l-F(s'))*' 


(1) 


where s = size (i.e., genomic content in bp) of DNA; R(s) = multi-track remaining-size distribution (Sachs et al, 
1999) for DNA fragments; I(s) = unity for s<S and zero otherwise (telomere cutoff (Laughton et al, 2004)); A = 
cluster intensity (Sachs et al, 1999) = number of one-track DSB clusters per unit size of DNA (A is proportional to 
the dose); M = average multiplicity (number of DSB) of an individual DSB cluster, induced by one track on one 
chromosome; and F(s) = cumulative distribution for one-track DNA fragment sizes. R(s), which completely char- 
acterizes the multi-track fragment size patterns (Sachs et al, 1999), can be fit to a high-dose PFGE experiment. 

Thus Eq. (1) expresses full information about the multi-track response in terms of the dose parameter A and of 
dose-independent, one-track quantities (M and F(s)). The proof of the equation is long but straightforward 
(Sachs et al, 1999). 

2.2. Monte-Carlo simulation of DNA damage 

To calculate the one-track quantities M and F(s), we introduce (except in the low-LET limit) specific models 
for the track and chromatin geometry. Our limit of resolution is set at 2 kbp for the smallest discernable feature of 
DNA (Figs. 1, 2) so that we can model the whole human genome. We consider the radial distribution of dose and the 
distribution of electron energies in the average track model (Cucinotta et al, 1999; Kobetich & Katz, 1969). In this 
model, the distribution in the number of electrons ejected, /?, with energy, co and solid angle, Q produced by ion 
interactions with target atoms (i.e., the double differential cross section for electron ejection) is combined with the 
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average transmission properties of electrons to obtain the spatial distribution of electron dose as a function of radial 
distance from the ion’s path. The radial dose due to ionization is then described by (Kobetich & Katz, 1969) as: 




(2) 


In Eq. (2), E is the residual energy of electrons after penetrating distance t, and rj(t,co) is the transmission probabil- 
ity that an electron with starting energy co penetrates to a depth t. Reference (Cucinotta et al, 1995) describes the input 
functions necessary to evaluate D in Eq. (2). Cross sections for electron production from protons are scaled to heavy 
ions using effective charge, which is expected to be accurate for high-energy ions. We follow (Brandt & Ritchie, 1974) 
for the dependence of the excitation term, D exc (t), on radial distance / in terms of the charge number Z and velocity with 
d(P)=(p/2)-hc/(2n(0 r ), where c is the speed of light,/! is the relativistic velocity factor, and h is Planck’s constant, and 
we use co r = 13 eV for water. The constant C,, xc is adjusted so that the model provides a good estimate of the LET 
(Cucinotta et al., 1999). 


exp (—t/ 2d 
2 

t 

This track structure model can be combined with a chromatin model to calculate the one-track action. In our 
NASARadiationTracklmage model, a random walk structure of chromatin is used as the basic DNA geometiy. We 
divide a chromosome into a large number of evenly spaced monomers, with each monomer being the same size, typ- 
ically « 2 kbp. The spatial geometry of the chromatin is approximated by a random walk on a cubic lattice. Each lattice 
site is a node on the cubic lattice representing physical space; a monomer is located on one such site. The lower limit of 
resolution in the model corresponds to one monomer size, below which the structure of chromatin is not specified and 
DSB distributions are not modeled. 

The track radial energy profile (Eq. (2)) is used to generate a stochastic process of DSB formation on human 
chromosomes, given by the probability *¥ of having at least one DSB on a monomer and given the monomer distance t 
from the track center in the XY plane ( i.e ., in the plane perpendicular to the direction of radiation). This probability is 
given by (Ponomarev et al, 2001a) as: 


D m {t) = C m {z,P) 



q/ = l_ e -e°W (4) 

where D(t) is the local dose (averaged over a small volume at the distance / from the track center). The quantity Q 
[Gy 1 ] is a track efficiency parameter; it is the only adjustable parameter in our model and is defined more specifically 
in Eq. (5) below. 

Despite large variations in D(t) for the high-LET radiation, we found that in a typical simulation of the 
PFGE experiment QD(t) « 1 in any location within the cell nucleus, even at the sites close to the track center. In 
this approximation ¥~QD(f), which corresponds to the fact that if the dose in any small physical volume inside the cell 
nucleus increases, the number of DSBs produced there increases proportionately. 

We can easily apply this formalism to the low-LET radiation by putting D(t) to const (and equal to the dose of the 
experiment), Eq. (4) in this case translates into a Poisson statistics, or random breakage (Cook & Mortimer, 1991). The 
stochastic process generated by Eq. (4) can thus be used for both low- and high-LET radiations. 

The algorithmic implementation of the RLC formalism can be expressed by the following two equations relating 
the experimental data and the formalism parameters: 


Q = DSB yield x DNA size / N 

(5) 

/l = QD/M 

(6) 


In Eqs. (5, 6), D is the dose of the experiment in Gy; DSBjyield is the number of resolvable DSBs per unit dose 
and unit size, obtained in a PFGE assay; DNA_size is 140 Mbp, an average size for a human chromosome; and N is the 
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number of monomers (=70,000) in an average chromosome. Another meaning of our adjustable parameter Q is the 
average number of resolvable DSBs per monomer per unit dose; that is, a local DSB yield. All radiochemistry involved 
in the creation of DSBs by ionizing radiation is encapsulated into the quantity Q, which will, of course, depend on 
track radiation efficiency but should not depend on dose (Eq. (4)). Eq. (6) shows how the dose determines k, the 
average number of DSB clusters per chromosome, or the cluster intensity, which is also the number of tracks that 
interact with a chromosome, each track having produced at least one DSB. Mis the average cluster multiplicity, or the 
average number of resolvable DSBs created on a chromosome in the one-track action (Sachs et al., 1999; Ponomarev & 
Sachs, 1991). The RLC formalism relates the multi-track fragment-size distribution function that can be obtained an- 
alytically from Eq. (1) and the one-track action (Sachs et al . , 1999) that is simulated by our algorithm using Eqs. (4-6). 
Our model is validated by fitting the DNA fragment frequencies obtained with our algorithmic and theoretical approach 
to the experimental data. 

Another justification of our approach to the prediction of DNA fragment sizes and the DSB yield is that the 
DSB yield is never directly measured in a PFGE assay: some approximations are used to infer its value (Lobrich et 
al., 1996; Hoglund et al., 2000; Radulescu et al., 2004). Our model provides a systematic approach instead of im- 
plicit extrapolations that PFGE experimentalists do for small fragments and the DSB yield inferred from the number 
of fragments. We predict the DSB yield from Eq. (5), once the optimal value for Q is found by / -fitting the model to 
the experimental data. The reason our DSB yield is more systematic is because it is based on a more detailed analy- 
sis of the track structure and the DNA geometry. We presented the /"-fitted Q values for a variety of experimental 
situations elsewhere (Ponomarev et al, 2001a; Ponomarev et al, 2001b). These values of Q can be used to predict 
the DNA fragmentation and DSB yields at any dose for any ion presented in Table I. The low dose extrapolation 
can be made with our model, and the DSB yield and DNA fragmentation can be predicted for a Microbeam 
experiment with a fixed number of particles per cell. 


Table I. The intensity parameter Q of the stochastic process of DSB formation. The estimated 
values for the intensity of the stochastic process of DSB formation Q(Eq. 1) are shown in a variety 
of experimental situations. Data are grouped into dose intervals and are taken from different groups 
(based on the data from (Lobrich et al., 1996)*; other entries based on the data from (Hoglund et al., 
2000); newer data* are based on a model fit to (Radulescu et al., 2004)). Some discrepancy between 
experimental groups was observed, e.g., Q=2.46*10 5 Gy -1 for (Lobrich et al., 1996), and Q=1.7*10 5 
Gy -1 for (Hoglund et al., 2000) for the same radiation (N ions at 80-97 keV/pm), and the same dose. 
Q is dose-independent and near LET-independent. 


QxlO 5 

Fe, 150 
keV/qm 

N, 80-97 
keV/qm 

N, 125 
keV/qm 

N, 175 
keV/qm 

N, 225 
keV/qm 

He, 40 
keV/qm 

60 Co- 

photons 

<0.5 

keV/qm 

150- 
200 Gy 

3.71±0.01* 

2.46*;1.70 

±0.01 

1.81±0.01 

1.43±0.01 

1.41±0.01 

1.83±0.01 

1.2±0.01 

100 Gy 


1.72±0.01 

1.78±0.01 

1.56±0.01 

1.46±0.01 

1.67±0.01 

1.2±0.01 

50-80 

Gy 

3.87±0.01* 

2.49*;1.64 

±0.01 

1.78; 1.81 
±0.01 

1.44±0.01 

1.47±0.01 


1.08; 

1.25;1.19+ 

±0.01 

30 Gy 



1.78±0.01 




1.25±0.01 




2.3. Model of large DNA loops 

Chromosomal loops arrange chromosomes into micelle-like structures, or the DNA “factories,” where gene 
transcription takes place (Miinkel et al., 1999; Marko & Siggia, 1997). The loops we simulate are different from the 
giant loops on Mbp scale hypothesized in the multi-loop-subcompartment model of DNA loops on the Mbp-Gbp 
scale (MLS) model (Miinkel et al., 1999). We do not model smaller loops (~ 20 kbp) either, because they are 
below the lower limit of resolution of the PFGE experiment. The lower limit of the resolution in our model is the 
Kuhn length (de Gennes, 1979), which is based on the statistical properties of chromatin from 20 kbp to the 1 Mbp 
range. Since we determined this length as 2 kbp, it is, therefore, the size of a monomer in our RW. 

We model the proposed loops on the scale of 120 kbp (Miinkel et al., 1999). In our model, each chromosome 
begins with a random stretch of 120 kbp, then 10 consecutive loops of the size 120 kbp are simulated, then goes the 
next stretch of 120 kbp of DNA and so on (Fig. 3b). It has been shown that the end-to-end distance between two loci 
(monomers) on DNA, separated by £ kbp, scales like an RW within a certain genomic size range. In (Miinkel et al., 
1999) it is argued that there are two linear regimes for the end-to-end distance scaling, which can be explained away 
by assuming that under -Mbp genomic sizes we have an underlying chromatin RW, and for distances above -Mbp 
scale we have an RW of larger aggregates of chromatin, which are suspected to be “giant” loops (Sachs et al., 1995; 
Miinkel et al., 1999; Marko & Siggia, 1997). In this work, we rely on data in Ref. (Miinkel et al., 1999) and assume 
that the loops are on average 120 kbp in length, and that there are about 10 such loops in a bundle called a polymer 
micelle (or rosette, Fig. 3b), with each micelle being separated by a run of chromatin of 120 kbp. 

Our code simulates each loop by a brute force approach, i.e., no artificial rule to close a RW to form a loop is 
enforced. Each loop has 60 monomers. The first loop monomer is placed at a random (x,y,z) point inside the cell 
nucleus. The program creates an RW shooting in an arbitrary direction (Fig. 2). Then the RW evolves for 60 steps 
subject to volume restrictions imposed by the nucleus membrane and the chromosomal domains. If on the last step 
the RW does not return to the origin (x,y,z), it is discarded. Only when an RW is generated that closes onto itself can 
we say that it is a statistically unbiased loop and a thermally equilibrated one (Kremer & Grest, 1995). A collection 
of such loops is in thermal equilibrium, because no artificial restrictions were used to close RW onto themselves, 
besides the restriction of the nucleus membrane and the domain boundaries. Such simulation is very central pro- 
cessing unit (CPU) intensive, as a significant fraction of computational time is wasted to form a single loop. 

2.4. Domain generation 

The generated domains are spherical with some overlap allowed. Their volumes are proportional to the sizes 
of chromosomes. Each domain is checked for intersection with other domains, and this statistics is recorded. The 
x,y,z location of the center of a domain is simulated and placed randomly within the nucleus. Flowever, domain 
centers cannot be too close to each other to provide for only partial overlap between domains. The diameter of each 
domain is calculated proportionately to the size of the corresponding chromosome. For instance, the domain for the 
first chromosome has the radius of 3.9 pm, while for chromosomes #46 the radius is only 2.3 pm. Forty-three 
percent of the nucleus volume is occupied by the domains. The domain center locations can change during 
the Monte-Carlo simulation, but the diameters remain the same. 

2.5. The implementation of the numerical model 

The DNAbreak model is now available as a part of the interface (NASARadiationTracklmage, copyright 
by Universities for Space Research Association (USRA)) written in Visual C++ 6.0, Visual FORTRAN 6.0, and 
OpenGL (open graphics library). This GUI is designed for experimentalists who will want to apply the model to 
their experimental data. All input data corresponding to experiments can be entered through dialogue boxes. The 
program used for DNA simulation, DNAbreak, was combined with Visual FORTRAN code previously employed 
by (Cucinotta et al., 1995 and 1999). These tools are now combined into one master interface. Data from other 
models such as chromosomes domains and stochastic track structure (Nikjoo et al., 2002a) can be input from data 
files. The GUI has three-dimensional (3D) displays of chromosomal domains, and will allow users to set up sim- 
ulations through a single dialogue box that contains all parameters and options needed to set up the track properties, 
DNA and nucleus geometries, and experimental conditions. The track, the nucleus, the chromosomal domains, 

DNA loops, DSB (Fig. 3c), and chromosomes can be visualized through an OpenGL graphics window. This 
interface will help researchers to plan future experiments, and will be a reliable tool in predicting the out- 
come of envisioned experiments. 
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Figure 4. Comparison of the model to the data. The 
fragment-size distribution function (% DNA) is plotted vs. 
the fragment size in Mbp. Circles are experimental data 
(Hoglund et al., 2000; human cells, N ions, dose=50 Gy, 
LET=225 keV/pm). The NASARadiationTracklmage model 
calculation is shown as a solid line. It fits data better than 
the random-breakage model (dotted line). 
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Figure 5. Comparison of the model to the data. The 
fragment-size distribution function (% DNA) is plotted 
vs. the fragment size in Mbp. Circles are PFGE data 
(Lobrich et al., 1996; human cells, Fe ions, dose=189 
Gy, LET=151 keV/pm). The least-square fit obtained 
(the solid line) determines Q (=4.94x1 O^Gy -1 ). The 
corresponding DSB cluster density is A=0.224 Mbp -1 . 
The Fe ions are the most significant component of 
galactic cosmic rays (GCR), and will contribute the 
majority of DNA damage. 


3. Results 


3.1. Prediction of the fragment-size distribution function 

Our model predicts the mass distribution of broken DNA fragments (Fig. 4). The graph shows, on a log-log 
scale, the mass fraction of a DNA pieces within a given size range vs. the average length of a fragment within this 
range. High-LET particles produce more relatively small DNA fragments, which is indicative of more concentrated 
and, therefore, less reparable 


biological damage. For instance, X 
rays would have produced fewer 
fragments on the left side of the 
graph. Flere we show the good 
correspondence of the model data to 
the experiment, the quality of 
radiation to produce DNA fragments, 
Q ( 1/Gy), the average density of 
fragments per genomic length 1 
(1/Mbp), the yield of the DSBs, 
which are responsible for DNA 
fragmentation in # DSBs/Mbp/Gy 
units. 

In Figs. 4 and 5 
we give an example of 
N AS ARadiationT racklmage 
prediction of the fragment-size 
distribution function based on the/ 2 - 
fitting of the experimental data. The 
model parameters are given in the 


Table II. The estimated values for the DSB yields, 
#DSB/Mbp/Gy, are shown in a variety of experimental situations. 
Data from (Lobrich etal., 1996( (subscript 3 ), and (Hoglund etal., 
2000) (no subscript) were used. Optimizing DSB yields does not 
lead to very different values, indicating consistency between 
simulated data (DSB yields and the fragment-size distribution 
function) and experimental data. NASARadiationTracklmage- 
estimated DSB yields are recommended for high-LET radiation 
as more mechanistically justified. 



figure captions. For comparison. 
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results of the random-breakage model are also given: the random-breakage model fails to give a useful prediction for 
the fragment-size distribution function of high-LET radiation. The discrepancy between the two distributions is less 
pronounced as LET decreases. The DSB yields for other simulations are summarized in Table II. 

The yields did not differ by more than 5% from the previously estimated ones. However, we suggest the 
NASARadiationTracklmage model as a robust tool to estimate these yields, as the model was shown to produce 
high-quality fits for various LETs. 


3.2. The dependence of the track efficiency parameter Q on dose and LET 


The values of Q for a number of experimental situations 
of the RLC formalism Q should be independent of dose. 
Intuitively, this can be seen from the fact that the DSB 
production efficiency of one track should not depend on 
whether other tracks are present. This property is con- 
firmed from the numerical fits to the PFGE data (Figs. 4, 

5). 

3.3. Simulation of the Microbeam experiment 
with a fixed number of particles per cell 

In this work our model was updated to include the 
simulation of an experiment with a fixed number of 
particles per cell, which we call here the simulation of 
the Microbeam experiment (Randers-Pehrson et al, 

2001). We intended to give a model representation of 
DNA breakage as it would occur in the Microbeam 
experiment. In the simulation, we choose a certain 
number of particles n hitting a cell with a randomly 
distributed impact parameter/), which originates from 
the imprecision of the particle collimator (Randers- 
Pehrson et al, 2001). The focus of the beam is char- 
acterized by the parameter a, which is the width at half 
maximum for the distribution of the impact parameter/). 
The number of tracks n can be any in our code, and in 
this work we demonstrate the results for situations with 
1, 2, and 16 tracks per cell, which correspond to real 
experiments (Randers-Pehrson et al, 2001; Sawant et 
al, 2001). We used the Q values obtained from 
simulation of PFGE experiments (in Ponomarev et al 
2001b, and in Figs. 4, 5) with a homogenous radiation 
field and a fixed dose. The number of particles in those 
was Pomon-distributed in the plane perpendicular to the 
incident particles direction. Simulation of PFGE 
experiments calibrates Q , which allows us to predict 
fragment sizes for other situations, such as the case when 
one has a fixed number of particles per cell with a non- 
homogenous distribution of the impact parameter 
characterized by a. We analyzed the theoretical data for 
their dependence on the most important parameters 
characterizing the quality of radiation and the geometry 
of DNA configuration in cell nucleus as well as the 
shape of the nucleus. 

Fig. 6 shows that LET is the most significant factor 
determining the shape of the number of DSBs per cell 
per track distribution, the shape of the nucleus is second, 


are given in (Ponomarev et al., 2001b). By the construction 



0 5 10 15 20 25 

n DSB per track 


Figure 6. A track striking a cell, as it happens, for 
example, in a Microbeam experiment, can produce any 
number of DSBs with different probabilities. Too many 
DSBs are highly unlikely; on the other hand, no DSBs 
can have a sizable probability for not-so-high-LET 
particles such as 40-keV/pm a-particles. Other geometric 
factors influence this probability, as we can see from this 
figure. Each particle has its own most likely number of 
DSBs and its own average number of DSBs. This figure 
shows that high-LET 225 keV/pm N ions are more effi- 
cient in the production of DSBs (hollow squares), and that 
LET is the most important factor influencing the number 
of DSBs per track. The group of curves above the hollow 
cubes is the action of 40 keV/pm He particles. Hollow 
circle: flat nucleus (2:2: 1 /2 main axes ratio), defocused 
beam (o=2.5 pm); filled triangles: round nucleus, defoc- 
used beam (o=0.17 pm); hollow triangles: round nucleus, 
a=2.5 pm, chromosomal domains imposed; filled rhombi: 
round nucleus, o= 0.17 pm; filled boxes: round nucleus, 
a=0.17 pm, domains imposed. The curves show impact 
on the nucleus geometric parameters and other param- 
eters on the efficiency of DSB formation. Flat nucleus hit 
by a moderate-LET a-particle has about 25% of “misses,” 
or no DSB per track per cell (empty circles). A round 
nucleus and a focused beam setup result in a Poisson- 
like distribution of the number of DSBs with the maximum 
probability for 2-4 DSBs per track. The effect of domains 
is not detectable. The least important factor after the 
beam focus and the flatness of the nucleus is the effect 
of chromosomal domains. 


the impact parameter (the focus of the beam) third, and the presence of chromosomal domain is the least significant 


factor. We came to a conclusion that DSB production per cell per Gy is the same for a constant dose by different 
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LET. This agrees with the established fact that the DSB radio-biological efficiency (RBE) is almost LET- 
independent. 

However, depending on LET, a single track can have a totally different impact on a cell. The average dose, 
based on the energy imparted to the nucleus volume, and averaged over that volume, depends on the type of inci- 
dent particles. If one averages the energy imparted over the cell volume, one can calculate the action of an N ion 
particle at 255 keV/pm corresponds to a dose of 0.43 Gy (this is similar to the estimates of the dose delivered by 
a single Fe ion mentioned in the Introduction), and a 40-keV/pm a-particle corresponds to a dose of 0.08 Gy, both 
for o=0.17 pm. A dose associated with one particle corresponds to the dose, as it is reported in an experiment with 
a random field of particles, if each cell gets one particle on average and the result is averaged over many cells and 
over cell volumes. Today technology allows irradiating a vast amount of cells, each in a similar manner (Randers- 
Pehrson et al, 2001). The reason why the distribution of the number of induced DSBs in Fig. 6 is so much wider 
for 225-keV/pm N ions is that there are many more DSBs on average per particle per cell for this particular ion. 

The DSB yield per Mbp per Gy is about the same for all ions. 

Fig. 6 shows that the particle LET has a significant impact on the distribution of the number of DSBs per 
particle per cell in a simulated Microbeam experiment. An important finding is that, for not so high LET, a particle 
hitting a cell does not necessarily produce a DSB. This is true, for example, with the He ions (the first point on each 
curve except the curve for N ions). We believe that this finding is important for the interpretation of the experiments, 
as we observe that, in 7% of all cells, no DNA breakage can be observed even if there is a direct hit by a 40-keV/pm 
a-particle. Higher-LET N ions lead to a practically nonexistent “no DSB situation” and a higher probability to have 
several DSBs per particle per cell. It has been asserted that heavy ions are not really “high-LET” radiation because 
of their wide penumbras, which sets the ratios of certain chromosome aberrations (Sachs et al, 1997) to be different 
from high-LET predictions for this biological endpoint. We see from our numerical experiments that a heavy ion is 
more of a combination of high-LET radiation focused closer to the track core, and low-LET radiation in the penum- 
bra that can extend beyond the nucleus. The overall effect is, however, a non-Poisson fragment-size distribution 
with higher clustering of DSBs within nucleus (Euclidian space) and along chromosomes (genomic dimension), 
as the radiation at the track center has the most significant effect. 


3.4. Analysis of the impact of the 
chromosomal domains on DNA 
fragmentation 

The simulated domains do not seem to 
have a large impact on the DSB distribution; 
this can be seen from Fig. 6. In Microbeam 
experiment simulations, we used our model to 
deduce data on the DSB number distribution 
that would be difficult to obtain experimentally 
because of the small number of DSBs produced. 
If one particle, such as a 225-keV/pm N ion or 
a 40-keV/pm a-particle, hits a nucleus, then 0, 

1, or more DSBs can be produced. The higher 
LET, the less likely 0 DSB (no DSB) will be 
created. For higher LET particles, we should 
see a peak in the DSB number distribution for 
larger numbers of DSBs produced per cell per 
particle. The average number of DSBs per 
particle could be low (about 4 DSB for 40- 
keV/pm a-particle). Notice, however, that the 
overall DSB yield in our numerical experiment 
corresponds to the ones reported in (Kremer & 
Grest, 1995), which are 7.8-8.2xl0~ 3 #DSB/ 
Gy/Mbp. This correspondence follows from our 
calibration of the model (Figs. 3, 4). Reported 
DSB yields for higher fluences are not defined 
the same way as the DSB yields per particle per 
cell. Our model calculates the yield of 
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Figure 7. Cumulative probability of fragments of length 
<t. The lengths shown correspond to only a fraction of all 
fragments. The DNA loops (dotted line, He 40 keV/pm) as 
seen from this figure will increase the frequency of shorter 
fragments in comparison to a simulation with no DNA loops 
(solid line). Higher LET (the dashed curve, 150-keV/pm Fe 
ions) leads to the same effect but to a smaller degree. This 
figure suggests that the effect of the DNA loops in combina- 
tion with high-LET ions, such as 40-keV/pm He ions, is the 
prominent excess of short fragments (less than 0.05 pm in 
length). The dot-dashed line is the data for the same He ions, 
but in the DNA domains imposed. The domains do not have a 
significant influence on the number of small fragments. 


#DSB/Gy/Mbp because we can calculate the average dose in the cell nucleus. For instance, one a-particle deposits a 
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small amount of energy on average corresponding to a dose of only about 0.08 Gy for a focused beam. The 
combination of the small dose and the small number of DSBs leads to a DSB yield (#DSB/Gy/Mbp) that is similar 
to the one reported in the experiment. 


3.5. Geometric parameters controlling DNA breakage 

The important conclusion that we can derive from Fig. 6 is that the DSB number distribution shows drastic 
change for higher LET, but for the same LET shows much smaller sensitivity to nucleus and DNA geometry. The 
curves for 40-keV/pm a-particles show that a 


flat nucleus would have a stronger impact on 
this distribution. The flat nucleus is defined as 
a disk with the smallest axis parallel to the 
direction of radiation (the principal axis ratio of 
the nucleus given 

by ax 2 + by 2 + cz 2 =(6 jum) 2 , 



cannot make any conclusion on the degree of 

the influence that the presence of 

the chromosomal domains, or the change in the 



impact parameter <r, should have, because those 
effects are small and cannot be discerned even 


# DSB 


for high Monte-Carlo sampling (-20,000 sam- 
ples). However, this is in an intuitive agreement 
because the chromosomal domains should only 
affect the locations of telomeres, which does 
not affect much the fragment-size distribution 
(Sachs et al., 1999; Ponomarev & Cucinotta, 
2006). 

3.6. The effect of loops 

RW alone cannot account for the 
complexity of DNA organization in the cell 
nucleus. DNA loops have a strong impact on 
the number of small fragments (Fig. 7). The 
DNA loops lead to a preponderance of smaller 
fragments (Ponomarev et al., 2006). This is in 
agreement with experimental work that points 
out that the high-order organization of DNA 
is accountable for the excess of the observed 



smaller DNA fragments (Radulescu et al., 
2004). 

3.7. Justification of the penumbra 
cutoff 

A numerical datum, in which the cut- 
off parameter for penumbra was varied with 
all other parameters constant, shows near- 
independence of the simulated track width. For 
N ions, the extent of the penumbra can be quite 
significant. The majority of numerical experi- 
ments were performed at the cutoffs 5 pm, 
which is more than sufficient for narrow tracks 


Figure 8. The simulated extent of the track penumbra. To 
make a single track simulation faster, a smaller cutoff on the 
extent of the track penumbra is beneficial. The cutoffs of 3 pm 
and 15 pm for 225-keV/pm N ions do not change significantly 
the cumulative DSB number distribution for one track per cell. 
The solid squares correspond to the cutoff of 3 pm and the 
hollow squares to the 15-pm cutoff (Panel A). Panel B is the 
data for the same ion and the same penumbra cutoffs (empty 
squares are 3 pm, and solid squares are 15 pm). Panel B 
shows the cumulative probability of a DSB to be anywhere at 
distance t not exceeding distance r in pm. Our choice of the 
cutoff parameters did not affect this probability. The cutoff of 
5 pm was taken for other simulations as the optimal one. 


of a-particles. Very large penumbra cutoffs will 

slow down the simulation, as we need to simulate the physical cylindrical space of the radius Rp+Rc (penumbra 
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plus cell nucleus). The sufficient extent of the 
penumbra cutoff was determined from a series of 
numerical experiments. We tested various cutoffs 
for N ions as their penumbra can extend up to 15 
pm. However, for cutoffs of 3 and 15 pm, we 
have not noticed any systematic deviation for the 
number distribution of DSB in one track per cell 
experiment (Fig. 8). For all other data, the cutoff 
was 5 pm. The deviation between the simulated 
curves did not exceed the error of the numerical 
experiment. This is agreeable with the picture of 
drastic drop-off of energy distribution in ion 
tracks (Cucinotta et al. 1995 and 1999). 



t, nm 


3.8. The effect of the track structure 

The cumulative probability to have a DSB 
is shown as a function of a distance <t from the 
track center (Fig. 9). This function is defined as 
a probability for a given DSB (regardless of the 
total number of DSBs) to be at a distance <t. 

Three cases are considered (1-3). DSB induction 
from 40-keV/pm He ions with (case 1) and with- 
out (case 2) DNA loops have same spatial pattern 
(two upper curves) of DSB formation. The DNA 
geometry has no impact on the DSB distribution 
in the Euclidian space (assuming DNA density is 
homogenous, but not necessarily constant, within 
the nucleus); only the track geometry does. DNA 
loops only influence fragment statistics, not the 
DSB location in the Euclidian space. The track 
structure for the 150-keV/pm Fe ions (case 3) is 
quite different from that of the He ions. DSB can 
form further from the track center for Fe ions than 
for He ions. The overall DNA damage from Fe 
ions is more spread out in the cell nucleus, even 
though the total number of DSBs per cell per 
track is higher (Table II). We found that He ions 
induce fewer DSBs but are more concentrated in 
the Euclidian space, while Fe ions induce more 
DSB, but some of them are more distant from the 
main grouping of DSBs. 


Figure 9. The cumulative probability to have a DSB at a 
distance <t from the track center. This function is defined as 
a probability for one DSB, if such to occur from an ion, to be 
at a distance <t. Three cases are considered (1-3). DSB 
induction from 40-keV/pm He ions with (case 1) and without 
(case 2) DNA loop leads to the same spatial pattern (two 
upper curves). The track structure for 150-keV/pm Fe ions 
(case 3, no loops)) is quite different from He ions, as some 
DSBs can form further from the track center for Fe ions. 
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3.9. Q (the intensity parameter of the 
stochastic process of DSB formation) 
is independent of dose 

Fig. 10 shows the values for the adjustable 
parameter Q as a function of dose. The numerical 
fits to experiments produce values for Q. which 
are, by and large, independent of dose. This is in 
agreement with the Eq. (4) that sets up the model, 
where the average number of DSBs per monomer 
is QD. The dependence on dose is already factor- 
ed in through D, and Q should be dose-independ- 
ent by definition. Fig. 10 shows this dose-inde- 
pendence, but indicates a slight dependence of this 
DSB production per Gy than X rays. 


Figure 10. The intensity parameter of the DSB stochas- 
tic process Q vs. the dose of the PFGE experiment. This 
parameter is the only model’s adjustable parameter defined in 
Eq. (4). The data are taken from Table I with the exact doses. 
This parameter determines V, the probability per Gy to have 
a DSB at a monomer. In this work Q is determined by x 2 -fit- 
ting of the model to PFGE data, and is dose-independent as 
seen from this figure. For low-LET 60 Co photons this param- 
eter is about 1.2 Gy -1 (open circles), and slightly higher for 
high-LET N ions (125 keV/pm, filled circles). 


parameter on LET. High-LET particles are more efficient in 
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4. Conclusion 

Our numerical model is now enhanced with chromosomal domains and DNA loops that are an essential 
component of the high-order DNA structure. The goal was to address the non-random distribution of 10-40 kbp 
fragments, for which these structures can play a role. We showed, however, that the chromosomal domains are not 
important for the specific biological endpoint of DNA fragmentation; but they certainly are important for other 
biological processes. 

In our model, the whole nucleus with the human genome is constructed with large scale details of DNA 
organization. Simulation at this scale reduces our ability to address the smaller scale. Our limit of resolution was 
set at 2 kbp, below which our polymer model of DNA cannot discern the details of DNA breakage. The scale that 
we have addressed, however, overlaps the fragment ranges produced in the PFGE experiment. 

Our model successfully predicts the distribution of larger DNA fragments and extrapolates data to lower 
doses and smaller fragment sizes. To account for a complex process of DSB formation in a small volume contain- 
ing a small piece of DNA, a single adjustable parameter, Q, was used to generate a point process of DSB along the 
DNA arc length C (Miinkel et al, 1999) with proper spatial and fragment-size distributions. The parameter Q is the 
intensity of the stochastic process of the DSB formation, and is the only adjustable parameter in the model that is 
determined by the model fits to the experimental data. We showed that the parameter Q only weakly depends on 
LET (Table I), which is consistent with the fact that DSB RBE is almost LET-independent. Our conclusion is that 
non-randomness in the DSB distribution resulting from high LET, and the clustering both in space and along DNA, 
is due to the superposition of DNA geometry and track geometry. This model takes into account two physical pro- 
cesses: a local process of DSB formation from a certain amount of incident energy in a small volume, as character- 
ized by LET-insensitive parameter Q; and a global process of many DSB and correlations among them, determined 
by DNA and track geometry on a large scale. 

The DSB yield in the model can either be determined from Q. or by a direct count of DSB in Monte-Carlo 
samples. Both numbers are very close to each other. The novelty of our approach is that the stochastic process 
generated with this parameter reproduces the clustering of DSB for high-LET and, therefore, provides a better 
suggested DSB yields than the ones proposed in (Stenerlow et al., 2003; Lobrich et al., 1996; Hoglund et al. 2000). 

With the help of our formalism, we extrapolated the fragment statistics to size 0 by estimating the cumulative 
probability to have a fragment shorter or equal to the monomer size (2 kbp). In the code, the genomic distance on 
the lattice between two DSBs corresponding to C<2 kbp is recorded as f(0), where/is the fragment-size distribution 
function. f(0) has a finite value corresponding to all fragments of the size l<2 kbp. This procedure is our worka- 
round for the problem of missing fragment sizes in the PFGE experiment below the minimal fragment size, and 
adds the DSB count for such small fragments to obtain the total DSB yield. More detailed data for fragments 
smaller than 2 kbp are obtained by other researchers (Nikjoo et al., 2002a). When the fragment-size distribution 
is determined for all sizes, the model predicts the DSB yield with higher precision than other models because the 
yield is based on a non-random DNA fragment distribution. As we noted, no significant discrepancy from the 
reported experimental DSB yields was found. We believe that this is due to the fact that high LET does not 
affect the DSB RBE, but it does affect the exact form of the fragment-size distribution. 

The problem with the PFGE experiment (Stenerlow et al., 2003; Lobrich et al., 1996; Hoglund et al., 2000) is 
that the data for small fragments (10-30 kbp) are absent due to limitations inherent in the experimental technique. 
The DSB yield can only be obtained from the fragment-size distribution function extrapolated for a wide range of 
sizes (from a few bp to ~Mbp). Experimentalists typically use the random breakage model for the extrapolation to 
smaller sizes. The random model, however, cannot account for high-LET effects in the case of HZE particles. Al- 
though our reported results did not show DSB yields noticeably different from the reported ones in experimental 
articles, the technique we propose is more consistent with the current knowledge of DNA organization and the 
track structure. 

We have come to the following conclusions for the effect of DNA loops and the chromosomal domains on 
fragment sizes: 

1) The chromosomal domains impose some topological restrictions on the location of telomeres within a cell 
nucleus and a limit on the mean square distance given between two monomers separated by a genomic distance C. 
The presence of the nucleus membrane and domain borders changes the scaling law with which the mean square 
distance depends on t. The square root of the mean square distance no longer linearly depends on f , as it would for 
an unrestricted RW. This is observed in the ~Mbp region of genomic distances (Sachs et al., 1995). This effect will 
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be vanishing for the small £’s (about 40 kbp and up to 1 Mbp) found with DNA fragments at high doses, especially 
from high-LET radiation. However, our conclusion is that the domains do not play any role in the fragment-size 
distribution function. 


2) DNA loops on 120-kbp scale lead to an even more pronounced excess of smaller fragments in the 10-40-kbp 
range than predicted previously by the DNAbreak model (Ponomarev et al., 2001a and 2001b). The presence of 
DNA loops does not affect the DSB yield. DNA loops should not affect the statistics of short fragments for low- 
LET radiation, such as X rays, as their homogenous radiation field leads to random breakage regardless how 
DNA is arranged. 

3) Other conclusions are: The spatial distribution of DSB is not affected by the DNA organization at all; it 
only depends on the track structure. DNA fragment sizes do depend both on the track structure and the DNA 
organization. A flat nucleus, as it often occurs when cells are plated, resulted in a less severe DNA damage per 
track per cell. A defocused beam in the Microbeam experiment results in a less severe DNA damage per track per 
cell, a-particles with moderate LET can result in cells with no DSB induced, even after a direct hit, with a 
significant probability. High-LET ions, such as N ions and Fe ions, would always result in some breakage. High- 
LET ions produce more DSBs per track per cell, but DSBs are less concentrated around a line in the Euclidian 
space because of wider penumbra, while a-particles result almost in a needle-like action. That is, if there are 
several DSBs induced by an a-particle, they are almost always located along a line in Euclidian space. 

The significance of our model is in its ability to predict the effects related not to the total number of DSB per 
Gy per cell but to their mutual locations in the nucleus. We showed how the spatial and genomic arrangements of 
DSBs are influenced by the loops and domains. The predicted effects will be important for the analysis of DNA 
repair and chromosomal aberrations. 
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